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i-G I Abstract. An original approach is suggested to analysis of full quantum Liouville 

flj ' equation for single electron (quantum particle) interacting with ideal phonon gas 

(harmonic boson thermostat). It is shown that under the thermodynamic limit this 
equation can be exactly reduced to a system of evolution equations connecting density 
C^ ■ matrix of the electron and its simplest irreducible correlations with amplitudes of one, 

C/3 , two, three and more phonon modes. Possible application of these new equations to 

explanation of the electron mobility 1/f-type fluctuations in semiconductors and other 
media is discussed. The special case of electron in static disorder is also considered. 
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(^ ■ 1. Introduction 

O I Here I leave for a time the problem of molecular Brownian motion and related mobility 

1/f-type fluctuations in classical fluids [ll[2ll3llllini[6l|71|8l|9] and turn to similar 
problems concerning random walks of electric charge carriers in conducting media and 

K^ ■ related electronic 1/f noise [TOl [HI [121 [13] • Although the second subject seems much 

H [ more interesting from practical point of view, the first one gave more possibilities to 

examine how and why the conventional kinetic approximations of statistical mechanics 
had lost so important phenomenon as 1/f noise. Now it is clear that the Boltzmann's 
conjectures are wrong, and in more correct approximations any particle of both dilute 
gas and dense fluid manifests unique random diffusivity and mobility (what, in essence, 
was predicted already in [H]). In the light of these achievements we may try now 
revision of kinetic models of quantum statistical mechanics, beginning, naturally, from 
the quasi-free electron in thermally vibrating crystal lattice. Or, in other words, electron 
in phonon bath (or, in generalized terminology, quantum particle in boson field). 
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Of course, real crystal lattices vibrate anharmonically. At that the anharmonicity 
gives rise to not only relaxation of any particular phonon mode but also 1/f-type 
fluctuations in rate of the relaxation [15]. These fluctuations can be transmitted to 
electron's relaxation due to phonons and thus to electron's mobility. 

But the electron-phonon interaction by itself causes relaxation of electron's velocity 
and energy. Hence, - as experience of investigation of the mentioned classical systems 
says, - it by itself can cause also 1/f-type fluctuations in the rate of this relaxation and 
consequently in diffusivity and mobility of the electron. Therefore it is reasonable for 
the beginning to consider electron in harmonic lattice, that is in ideal phonon gas. 

In this paper our purpose is to reveal a minimal set of relevant statistical 
characteristics which completely and exactly determine all effects of the electron-phonon 
interaction, and derive exact evolution equations for such characteristics. The results 
will present a necessary mathematical base for further investigations. 

2. The model and the problem 

Let us consider a single quantum particle interacting with a potemtial field formed by 
very many mutually non-interacting plane- wave quantized excitations. The latter will be 
called "phonons" while the particle "electron" , although it will be spinless, or "Brownian 
particle" (BP), since its interaction with phonons may enforce it to a random walk, at 
least in the thermodynamic limit and at non-zero temperature when the phonons can 
serve as thermostat with infinitely large total energy. 

2.1. Basic equations 

The Hamiltonian of the mentioned system will coincide with that of widely used simple 
models (see e.g. [Ml [HI [18]) of quasi-free (conduction) electrons in movable (vibrating) 
crystal lattice: 

H = He + Hr,h + -f^i 



ph ~r J^int i 



He = - — , Hph = 2_^ ^k O'kO'k , 

k 

Here Q is the system's volume, so that wave vectors in the sums are separated by 
reciprocal volume d^k = {2tt)^/Q , and 

c-k = 4 (2) 

for any wave vector k . 

Being interested first of all in statistics of the electron's, or BP's, walk, we must 
consider the von Neumann evolution equation, 

zhp = [H,p] = Hp~pH , (3) 
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for full density matrix of the system, p . If at that we are interested first of all in finite 
temperatures, T = /3~^ > , then it is reasonable to describe phonon dependence of the 
p in terms of coherent states and the remarkable "Sudarshan-Glauber P -representation" 

[iniEolEi] : 

P = /WHP(a.a-)<i„V„. (4) 

where a = {ak} , |«)(tt| = Ilfc I'^A;)(q^A:| , da* da = Ylk^^l^^k = Ylk^^^^dQak, 
and the "quasi-probability density" [El [20l [21] , P{a,a*), in respect to the electron's 
variables still acts as an operator. Taking into account that (for any k ) 

a \a){a\ = a \a){a\ , a"'' |a)(«| = {d + a*) \a){a\ , 
|a;)(a| a = (d* + a) \a){a\ , |a)(a| a^ = a* \a){a\ , (5) 

with 

d^^, d* ^ 



da da* 

one can transform EqlHl into 

thP=[H,,P] + J2 ^^ ^^l^l - "'^^'^) P + (6) 

k 

+ ^ E { 4«/.[e^'^ P] + Cka*^ [e-''\ P] + 49^Pe^^- - Cke'^^^OkP } 

Further, it will be convenient to treat an electron dependence of P in the coordinate 
representation, introducing function {r\P\r') instead of the operator P, with \r) 
denoting eigenstates of the electron's coordinate operator r . Then it is convenient also 
to introduce new independent variables: 

ly I ly 

X = ^— , Y = r-r', 

Ak = akt^"" , A*^ = ale-'''' (7) 

Designating the matrix element {r\P\r') by the same symbol P , 

{r\P\r') = P{t,X,Y,A,A*) , 
we obtain for it from Eqj6]the equivalent equation 



ihP = Vy 

m 



P + 



k 

+ Y,huJk{A*^dl-Akdk)P + (8) 

k 



with 



d „ ^ Q ^ A* ^ 

Ok = TT^ , Oi, 



X = 77T7 ) "^Y = t:;7T , ^k = 777- , <Jk 



dX ' ' dY ' " dAk " dAl 
The Fourier transformation J exp {—ipY / h) P d^Y / (2tt hY gives the Wigner represen- 
tation with P becoming literally "quasi-probability density" . 
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2.2. Initial conditions 

Being interested first of all in thermodynamically equilibrium electron's "Brownian 
motion" , it is natural to start from such initial density matrix, - e.g. at t = , - which 
says that the phonon subsystem has canonical equilibrium probability distribution with 
given temperature T = 1//3 while the electron is located in vicinity of some given 
point of space, - e.g. coordinate origin, - and has equilibrium Maxwellian momentum 
distribution corresponding to the same temperature. Such the density matrix can be 
expressed by formulas 

p{t = 0) oc \^){^\exp{-PHp^,) , (9) 

3/4 

\E'(r) = ( — -:- 1 exp 



f2mTY'^ ( mT .\ 



In order to rewrite this in the P-representation, notice [2T] that for any particular 
phonon mode 



(1 — e '^^'^) exp{—/3fluja^a) = / |«)(«| /{u, a, a*) da* da 



where 



f{uj,a,a*) = exp 

N{uj) ^ 



\a\' 



N{uj) 



(10) 



exp {f3huj) — 1 
Therefore initial condition (|9]) is equivalent to 

Pit = 0) = vl>(r)v]/*(r ') W ficok, a^, a*,) , (11) 

k 

which is normalized distribution. 

Obviously, such initial condition as ([9]) and fITT]) has serious defect: it neglects 
statistical correlations between electron and phonons inevitably arising from their 
interactions. However, this defect may be simultaneously its advantage: it helps to trace 
the development of the correlations during subsequent evolution directed by Eqs 13)16] or 
IHl Therefore, it is reasonable to accept such initial condition at least for a time. 

2.3. Marginal distributions 

The object of our direct interest is the marginal density matrix of electron, 

Po = n '^^'^ '^ ' (12) 

k 

- with Trfc being trace over states of phonon mode k , - or corresponding marginal 
distribution 



Po = I Pfl daldak = [ P Y[ dAldAk 



Electron in ideal phonon gas 5 

Alternative description of its evolution is presented by hierarchy of equations for the 
chain of various marginal, or particular, distributions 



Pn= f P H dAldA, 



k^Kn 

where Kn = {fci . . . fc„} is a set of wave vectors of n different phonon modes. 
Applying these integrations to Eqsj6] or |8l we come to hierarchy 

Pn = -VVxPn + 5Z [2fc + f^Ak + fiDk]Pn + (13) 

k e Kn 



+ /i ^ Ag Pn+i dA^dAg , 



where P„+i corresponds to -ft'n+i = {Kn, q} = {ki ■ ■ -kn^} , and we introduced 
operators 

V = Vy , 

m 

Lk = tiujk- kV) (AA - Aldl) , (14) 

Ak ^ {ih^r^ ( e'^""^ - e-^'^^/^) ( clAk - c^Al) , 
Du = {ih^^)-^e-'^'"\cldl-cA) , 



(15) 



and besides the small parameter 



with VLq being some fixed "microscopic" volume. 

It will be useful to treat these equations in terms of new functions Qn defined by 

-' n ^n I I Jk ) 

k&K„ 

with shortened notations 

h = f{uk,Ak,Al) = {7rNk)-'exp{-\Ak\yNk) , 

Nk = N{cok) = l/[exp{(3hwk)-l] 
They undergo hierarchy of equations 

Qn = -VWxQn + 5Z [-^-^ + l^^f" + l^fk^P>kfk]Qn + 

kdKn 

+ /i 5^ I jgAgQn+idAldAg = (16) 

= -VVxQn + 5^1 [2fc + l^Pk + l^Dk]Qn + 

keKn 



+ /" I] fqAgQn+idAldAg , 



where 



B, = Ak + (th^o)-' e-'"''/' ( CkAl - clAk )/N, (17) 
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Clearly, Qq = Pq , and initial conditions to EqsJT6] following from fITT]) are 

Qn{t = 0) = ^(r)^*(r') (18) 

Equations flT^ and flT^ are analogues of the BBGKY-like hierarchies of equations 
for BP in a classical ideal gas of atoms [H |2l [3l HJ [5] or generally in a simple fluid 
[21 in El El El 112] , with the last terms in ( 1T^ and (ITB]) playing role of "collision integrals" . 

The problem is again in investigation of properties of the BP's (electron's) 
distribution function Pq = PQ{t,X,Y) what follow from EqsJTSl or [TBI 

3. Thermodynamic limit and shortened hierarchy of evolution equations 

In order to grant to the BP (electron) large enough space for random walk and make 
the phonon subsystem to be as much as possible powerful thermostat, we must go to 
the thermodynamic limit fi — t- oo , yU — )■ . 

3.1. Preliminary discussion 

Notice that under this limit, firstly, perturbation of any particular phonon mode by its 
interaction with electron is vanishingly small, so that 

Pn^ PO n ^'^ ' Qn^ Qo = PO (19) 

k&K„ 

for all n . More precisely, the Eqs JSP^ and UHl along with initial conditions (llip and 
f lT8|) . do imply that 



f= E^'^'" IlA 



Ik 

.5 = J k 



where P*-**^ is a s -order polynomial of phonon variables Al, Af^ , with coefficients 
being some functions of t, X and Y . Correspondingly, 

oo 

Q^ = Y, l^' P^'^ , (20) 

s = 

where Pn is a s -order polynomial of Al, Ak with k E Kn and with coefficients 
representing some functions of t, X, Y and /i . 

At the same time, secondly, because of finiteness of limits 

- with dots replacing a "good" function of q (which is the case under natural restrictions 
on |cfcp), - summary perturbation of electron by all the phonon modes acquires at 
yU — )■ a constant value independent on /i and indifferent to exclusion (or neglect) of 
any finite set of phonon modes. Therefore, coefficients of all the polynomials Pn in 
fact become independent on /x . 



\3 
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These remarks show that, in agreement with flT^ . all P„ with n > give no 
essential information in addition to one already contained in Pq . Actually, at n = 1 
a portion of new additional information is presented by second term of expansion f l20p . 
i.e. by coefficients of the first-order polynomial P| . Indeed, according to the first 
of EqsJT3] and [16] (for Pq and Qo), at fi — )■ just P| directly and completely 
determines "collision integrals" in these equations (since /i^ ^ . . . -^ Qq f . . . d^q/{2-nY 
and, clearly, contributions from all Pi with s > 1 disappear in the thermodynamic 
limit). 

Similarly, any of Pn with n> 1 gives no new information in comparison with that 
already contained in coefficients of PI . Considering evolution of these coefficients, - 
with the help of second terms of flT3|) or flT6l) (at ra = 1 ), - one can see that at /i — >■ 
they are exactly determined, through corresponding components of collision integrals, 
by coefficients of polynomial Pg • O^; to be more precise, four coefficient of bilinear 
part of Pg . 

Continuation of such reasonings leads to conclusion that hierarchy of functions, 
which unambiguously determine evolution of the electron's distribution function Pq , 
in the thermodynamic limit reduces to coefficients of purely multi-linear parts of 
polynomials Pn (which is easy understandable: since density of phonon modes in A; - 
space tend to infinity, probability of strong excitation of one and the same mode becomes 
infinitely small as compared with probability of weak excitation of many different, let 
close in /c -space, modes). 

In order to confirm these statements, we have to extract from ( TT3l) or ( TT6l) a closed 
hierarchy of shortened equations for the mentioned coefficients. 

3.2. Derivation of shortened equations 

First, let us define them formally: 

1 



A„ = lim 



= A„Q„ , (21) 

A=A*=0 



n ^^^'^ ^« 

keK„ J 

where a^ can take one of two symbolic values, " + " or " - " , and we introduced new 
designations (more comfortable for further manipulations): 

^k ^^k , At ^Al , d-, ^ d, , dt ^ dl 

Applying such operation to {n + l)-th of EqsUB] (for Q„ ), consider different terms there 
separately. 

For the second right-hand term one readily obtains 

A„ LkQn = - icTk (wfc - kV) An , (22) 

where k is one of elements of the set Kn from f l?T]) . 

The third term needs in more careful consideration. Application of A„ to it gives 

AnAQn = -{th^/Q'o)-^akcl' X (23) 

X [e^i^y/^ - e-'^yi\\ + AT-i) ] K-xQn U,=a- =o , 
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where 



c, = c 



k 5 "-fc 



Cl =Ck 



and, clearly, A„_i corresponds to the set Kn-i = Kn Q k . 

Next, to express A„_iQ„ in ( l23l) through functions ( pTl) . we should take into 
account, firstly, that by the Qn-s definition 



/ Qn fq dA^dAg = Qn-l (K„ = K^^i © q) 



(24) 



Therefore 



^ n «rQ. 



A:e A"n-1 



/, c^A^c^A, = ^ n ^r Qn-l (25) 

Secondly, turning here all A^, At with k G -ft'n-i into zeros, in view of the expansion 
(12U]) for the square bracket on the left we can write 



<n-l K„_i ■ 



]^=^.=o + X]/x^p(^HA,^;) , 



(26) 



s = l 



where Ak„_i means any of Ak with k G Kn-i , while A any of all phonon variables 
with wave vectors from Kn , and p^'^\Ag, A*) is s -order homogeneous polynomial 
resulting from Pn^ ^ ■ 

Combining (1251) and (l26l) . in the limit ;U — ;■ we come to equality 

A„_iQ„|a,=a*=o = A„_i {Kn = Kn^i®q) (27) 

With its help EqES] yields 

KAQu =-{ih^or'crkcl' [e'''''/^-e-'''''/\l + N-')]An-i (28) 

(with Kn-i = KnQk). 

Now, consider action of A„ onto fourth right-hand term of (IT6ll : 

Ani^DkQn = lim /i2(i;iv^)"^e-*'=^/2 x 



M->.0 



/^ 



^ ( 49* - cA ) n ^."^ ^^ 



qeK„ 



A=A*=0 



Evidently, the square bracket here has a finite limit at /i — )■ , again thanks to ( 120|) . 
Therefore the total expression tends to zero, so that 

AnfiDkQn = (29) 

At last, consider action of A„ onto "collision integrals in f jTBjl . At that, we should 
use identities 

j UA,[...]dAldA, = J{-N,d;f,}A,[...]dA;dA, = 

= N,lf,d;[...]dA;dA, , 

j f,Al[...]dAldA, = N,jf,d,[...]dAldA, 



'■q ) 



(30) 
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where square bracket means arbitrary function of the phonon variables. Hence, colhsion 
integrals can be rewritten as follows, 



^n+l Qn+l = f^ / ^ / Jq ^q Qn+1 dA^clAg 



/i(2 



h^/Wo)-' Yl (e*^^/2-e-^«^/2)^^ 



X 



(31) 



qj^Kn 



and 



X / /g ( c*d* - Cgdg) Qn+i dA*gdAg , 

^n ^n+1 Qn+1 = 



X 



q^Kn 



X 



<7,e{+-} 



^q<^q 



l^ 



ir^;- n ^^^"^-1 



k£Kn 



Qj/\^Ciji\.Q 



Here y4x„ means any of A^ with k G -ft'n . According to fl20|) . the square bracket here 
behaves at /i — )■ in full analogy with that in ( 126|) . 



• • •J^K„=^|f„=0 



u=^*=o + $^/i^p^^n^.,^; 



(n+l+s) 



with polynomials p^'^\Aq^A*^ now produced from P^"^ 

Consequently, after performing the limit and transforming the sum over q into 
integral, we obtain 



ih 



J2 ^/c(e^''^/^-e-^'^^/2)iV,A„+i 



(2vr) 



(32) 



'^6{ + ,-} 

where Kn+i = Kn © q and a corresponds to q (so that one can replace dummy 
variable q by k^+i and its pair a by cr„+i). 

As the result of relations (^, (^, (^ and ([32D, we come under the 
thermodynamic limit to the following exact short-cut (nevertheless, infinite) hierarchy 
of evolution equations: 

A„ = - Wx A„, - z ^ o-fc {ujk - kV) A„, - 

keKn 

- {ih^/Wo)-' J2 ^k ^l' [e'"''^' - e-*'''^'(l + N^^) ] A„_i + (33) 






(T / c-'^(e^''^/2-e-^'^^/2)Ar^A„+i 



(2^ 



'^e{+,-} 

Here K^-i = Kn k and S„_i = S„ cr^ in the third term on the right and 
Kn+i = Kn q and S„+i = S„ o" in the last tern (colhsion integral), with S„ 
denoting sets of the two-fold indices paired with wave vectors. 
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Recall that all wave vectors in any of the sets Kn have appeared different one 
from another, so that "collision integral" in the last term, strictly speaking, does not 
include (infinitely small neighborhoods of) n points q G Kn . However, both the A„ 's 
definition (I2T1) and structure of Eqsl33] allow to extend functions A„ to coinciding wave 
vectors by continuity. 

Initial conditions to Eqsi33| what literally follow from (|9]) or f lTSj) . 

A„(t = 0) = (5„o^(r)^*(r') , (34) 

show that there is no need in any a priori assumptions about these functions. 

As for the electron's probability distribution function Aq = Pq , in essence, Eqsl33] 
altogether form an exact time- non-local kinetic equation for it. At that, its part of the 
initial conditions can be chosen arbitrarily. 

3.3. Stationary solution, equilibrium distribution, and thermodynamically improved 
initial conditions 

The Eqsl33] by their derivation describe an electron in infinitely large phonon (or 
generally boson) thermostat formed by continuum of modes each being initially in 
canonical equilibrium state. Therefore, firstly, stationary solution of EqsJ33| - to be 
denoted as A*^ , - must represent thermodynamically equilibrium distribution of the 
electron's momentum under fully uncertain electron's position, so that 

Vx At" = A:" = 



n n 



Such solution can not be strictly normalized to unit, - in respect to both momentum 
and coordinate of the electron, - but it can be normalized in respect to momentum only: 

Al'^iY = 0) = 1 

Secondly, time-dependent (strictly normalized) solution of Eqsj33] must tend to such 
stationary equilibrium solution in the sense that 

lim f And^X = A"^ 

t— >-oo J 

The functions A^^ contain complete information about equilibrium statistical 
correlations between electron's momentum and thermostat. Hence, if we want to take 
into account these inevitable correlations from the very beginning, we should replace 
initial conditions (j34l) by 

A„(t = 0) = Wo{X)Al" , (35) 

where a normalized probability density Wq{X) introduces initial localization of the 
electron in configurational space. Then time-dependent solution of Eqsl33]will highlight 
specific "historical" statistical correlations [21 H] arising between current state of the 
thermostat and summary path of BP (electron) during all the time interval (0, t) . It 
may be said that the "equilibrium correlations" describe momentary electron-phonons 
interaction while the "historical" ones describe a heritage of past interactions. 
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4. Equivalent representations and generating functional of electron-phonon 
correlations 

A suitable change of the variables A„ at n > helps to remove the formally arbitrary 
parameter Qq and fashion Eqs|33] to a more visual form. For instance, introducing, 
instead of A , functions D„ defined by 

-I— r C"^ 

An = Dn n Ska, , Ska = "T J^ r^ (36) 

k^Kr, nhJkNk^/ilo 

at n > and Dq = Aq , we transform Eqsj33] into 

Dn = - VVx Dn ~ lJ2 ^fc (^'^ - ^^) ^- - 

keK„ 

-I Y, Ok^k\e'^'"''Nk - e~'''''l\Nk^\)\Dn-x + (37) 



keKn 



l^gl / ^iaY/2 -iaY/2\ 7-1 "- ^ 



'^e{+,-} -^ 
with the same meaning of a and q as before. 



4.I. Kinetic equation and infinite hierarchy of irreducible correlations 

Physical and statistical interpretation of just made change becomes clear if notice, firstly, 
that, - as one can prove in the spirit of above consideration, - at different wave vectors 
definition (|2T|) is equivalent to 

^-= U Wk J^^o i^ /^" n ^r^ dAldAk (38) 

keKn " ^ -^ keK„ 

Secodly, if the electron was immovable ( "infinitely hard" ) then perturbation of all the 
phonon amplitudes would reduce to their constant shifts by 

— T'^ ^h n T 

^k = - ^ r^ = l^SkaNk 

riuJk V" 
Comparison of this expression with f l36p shows that 

Pn H ^ dAldAk , 

fcGA'n ^k 

that is Dn ( n > ) represent joint perturbation of n (different) phonon modes 
measured in natural units of their separate static perturbations. Under the Wigner 
representation of Dn 's electron dependence this relation can be rewritten as 

Dn{t,X,p,Kn,T.n) = D,{t,X,p) I \[ 4^) , (39) 

\keK„ ^k I x,p 

where the angle brackets mean conditional average of (quasi-classic) phonon amplitudes 
at given electron's position and momentum. 
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At that, our basic definition f l2T|) prompts that functions Dn {n > 0) indeed must 
be extended to (zero-measure subspaces of) coinciding wave vectors by continuity. 

From the Hamiltonian ([1]) it is seen that, in case of immovable ("static") electron all 
phonon modes would be statistically independent one on another, and all the statistical 
moments in (1391) factored into products of first-order moments (average values). In 
reality, however, electron's motion couples different phonon modes one with another 
making them statistically dependent. Therefore their statistical moments consist of 
not only average values but also second- and higher-order irreducible correlations 
(cumulants). They reflect both current and old interactions (flguratively speaking, 
"Correlations = Interactions + History"). The EqsJST] exactly describe influence 
of "current" onto dynamics and "old" onto statistics of electron motion (for detail 
explanations of nature of "historical" correlations see e.g. [71 [TOl [HI [121 [El [E]) • 



4-2. Generating functional 

It may be convenient to accumulate all the components of Eqsl37| into a single object, 
such as vector {Dq , Di , . . .} in the corresponding Fock space or the generating 
functional: 



oo _ „ n 

71=1 ' (T1...(T„ .7 = 1 



SO that 



Dr, 



n 



\ Sz„^{kj] 



V 



2 = 



In terms of this functional the whole inflnite hierarchy (ISTjl reduces to single equation 

6 






(40) 



where 



Lkc 
Bkc 

Akc 



VVxV + C{z, — \V 



icr [uk — kV) . 
\ck? 



-ikY/2 



{Nk + 1) 



(41) 



la- 



, ikY/2 



-ikY/2 



h^UJk{2Txf ' ' ' 

and V is electron's velocity operator deflned in ([T^ . 

Quite similarly we can introduce generating functional A for the correlation 
functions A„ and then write out similar evolution equation for it. Clearly, both they are 
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connected with V and EqJlO] through simple scale transformation of the test function 
z = z„{k) , 

A{^} = V{sz} , 

A = -\/VxA + £L^,^| A , (42) 

with above introduced multiplier s = Ska ■ 

The functional evolution equation ( l40l) is direct formal analogue of equations for 
BP in fluids investigated in [11 |2l |3l HJ |5]. Investigation of EqJlO] or EqJl2] may be 
subject of future works. At the rest of this one we confine ourselves by discussion of 
some consequences and modifications of the EqsJlU] and |121 

5. Generalizations and special cases 

5.1. Other particles and external potentials 

In presence of additional "external" particles, interacting with the phonon thermostat 
and distributed with given density z/(r) , the Hamiltonian ([1]) transforms as 

H ^ H + H, , Ha = -=^[c*i/*afc + CkVka\] , (43) 



k 

where 

ikr , 



z/fc = / e "'''iy{r)d-^r 

Then after repeating, for new Hamiltonian (1431) . the above expounded derivation of 
shortened equations, one easy comes to EqHO] with modified third right-side term: 

C ^ C' = C + ^ I d'^k z,{k) [iauk^k exp(mA;X)] , (44) 

with z/^ = i/fc , u'f^ = ul . This addition to evolution operator jC does not include the 
electron momentum related operators V and Y . Therefore it does not create additional 
electron-phonon correlations (instead producing momentum-independent shifts of mean 
values (Al) only). Hence, it can be removed by representing solution V to the 
modified evolution equation, 

V' = -VVxV + C'lz,-^\v' , (45) 



in the form 



V = V exp 



Y^ I z,{k)vle'''^^ d^k 



(46) 



Then a simple algebra turns the EqJlS] into 



i,=_,7v,I,+ ^Ml±Ml, + fL^^p, (47) 

in \ oz ^ 
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where r , r' = X ± Y/2 and 

is potential of effective (mediated by phonons) interaction of our electron with the 
"external" particles. 

This result is not significant in itself, for it follows already directly from the 
Hamiltonian (H3ll . with the help of boson operator amplitudes shifting 

Ak + ^^ ^ A, , AI + -^^ => Al 

Nevertheless, it is important for us since it verifies faultlessness of our derivation of 
EqsJ33]or EqsEZland thus EqsgDl 

5.2. Many phonon branches (boson thermostats) 

One of generalizations of the model ([T]), which may be necessary in its applications, 
is accounting for presence of several types, or "branches", of phonon modes, with 
different dispersion laws and/or couplings to the electron. Obviously, to perform such 
a generalization, it is sufficient to impart an additional index, - let "b" , - to the phonon 
amplitudes Ak,Al, frequencies Uk and coupling constants c^ , the sets of n such indices 
to the correlation functions A„ and Z)„ and, correspondingly, the same additional index 
to the test function z„{k) and all operators f l4T]) . so that 

CT b (J 

in EqsJlU] and other functional evolution equations. 

Equivalently, we can divide the corresponding total functional evolution operator, 
C , into sum of operators C'^^^ each relating to one of the phonon branches, so that the 
EqJlO] takes form 

V = -VVxV + Y,'^^'Azb,jA'^ (49) 

Advantage of this representation is in that it naturally extends to the case of different 
temperatures of various phonon branches (in other words, boson thermostats), which is 
obviously allowed by our above analysis. 

5.3. Application to electron in static disorder 

Let us consider a quantum particle ("electron") in Gaussian static random potential 
$(r) with zero mean value and correlation function 

($(ri)$(r2)) = [ \ck\^e'''^''-''U^k/{27i)^ (50) 



We assume that the electron's density matrix, pe , at some time moment, t = , is 
known, pe{t = 0) = peo , and want to know its mean value at t 7^ , 

nt) = (e-^^^/'^Peoe^^*/^') , (51) 
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where H = p'^/2m + $(r) is Hamiltonian of the system. 

In calculations of this quantity the so-called operator representation of random 
processes and fields may be useful [22]. For instance, any statistical moment of Gaussian 
field with auto-correlation ( l50l) can be represented as 

illHr,)) = (0|n$M|0) , 



where 

$(r) = -, 



is quantum boson field operator with coefficients satisfying ([2]), f2 — )• oo , and |0) is 
the "vacuum state" . Therefore instead of (1511) one can write 

now with H = p^ /2m + <l>(r) . 
From here it follows that 

mU-t) m = Tie Peo{Tr, e-^^*/^|0)|^)(^|(0|e^^*/n , (52) 

where \^) is arbitrary electron state and Trg , Tr{, mean traces over electron and boson 
degrees of freedom, respectively. Obviously, the expression in the braces is particular and 
peculiar case of the expression ( !T2|) arising from it under limits T — > and cj^ — )■ . 
At that, the first of these limits should be performed before the second, since transition 
to vacuum initial state of boson subsystem presumes that A^^ — > . 

To realize such ordered limits in EqsJ33| we firstly have to rescale functions A,j by 

^n = J^n II SkcTk J Ska = "ITTT /^ t 

instead of f l5B]) . then set A^^ = and after that w^ = . This results in equations 
Dn{t, X, F, K, S) = [- Vx + tY,akk]VDn{t,X, F, K, S) + 

+ i ^ afc e-""^'^ A.-1 (t, X, F, ir e fc, S e fffc) + (53) 

k&K 

(with designations slightly different from that in Eqsl37] but clearly equivalent to 
them). According to ^^ and (IS21) . if Z^„(t) is solution of these equations under initial 
conditions 

Dn{t = Q) = <5„o|^)(^| , (54) 

then 

(*|p,(-t)|^) = TrePeoI^o(t) (55) 

In particular, at peo = |^o)(^o| we obtain expression, 

(K^le^^^/'^l^o)!') = {^o\D,{t)\m,) , (56) 
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for mean probability of electron's transition between two given states \E'o and \1' . 

Because of the "Anderson localization" of a part of electron states in random 
potential [231 1211 125], the quantity (!56|l has non-zero limit at t — )■ ±oo depending 
on both \E'o and \E' . This fact means that stationary solution of Eqsj53] is not unique 
but strongly degenerated (consequence of that now the "thermostat" has zero internal 
energy and thus no ability to "thermalize" electron). These multiple stationary solutions 
can be used to represent such ("eigen") electron states which are correlated with the 
disorder (in contrast to a priory given states \E'o and \E'). 

5.4- Naive "weak- coupling limit" and conventional kinetic equation 

Returning to the case when T ^ and Uk y^ , let us assume that some sufficient 
smallness of the coupling constants |cfcp , e.g. in the sense of 

^ = J ¥^, J2^ '^ 'h ' 

gives us rights to truncate the hierarchy fl33|) or, equivalently, ( 137]) at its second level by 
setting A2 = or D2 = , i.e. neglecting third-order correlations (between electron 
and two phonon modes) and higher-order ones. Then, after transition to the Wigner 
representation, so that 

exp {±ikY/2) F{Y) =^ F{p t hk/2) , 

we come to equation 

D,{t,X,p) = -VVxD,{t,X,p) + J2j J^-j^J^ ^^ ^ 

hk"^ 

X {exp [-ia{ujk - kV + - — ) r ] x 

hk 
X [Nk Do{t-T,X-T{V-—),p- hk) - 

2m 

hk 

-{Nk + l)Do{t-T,X-r{V-—),p)] + 

2m 

hk"^ 
+ exp [-iaiuk - kV - 7;— ) r] x 

2m 

hk 

X [{Nk + 1) /^o(t -T,x-T{y + —),p + hk)- 

2m, 

hk 
- NkDoit - T,X - t{V + —),p)]} , (57) 

where V = p/m,. Further standard assumptions about spatial dependence (i.e. X- 
dependence) of the density matrix and standard reasonings (or "conjurations"), along 
with formal identity 

lim y / exp [cr ( . . . ) r] c?r = 2n5{...) , 
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help to transform f l57|) into what is called "kinetic equation", 



Do{t,X,p) = -VVxDoit,X,p) + (5^ 

+ TT / \ckm6{uJk~kV + 



h? J 2m 

x[NkDo{t,X,p-hk)- {Nk + l)Do{t,X,p)] + 

+ 5{ujk -kV- -— X 
2m 

d?k 
x[{Nk + l)D,{t,X,p + hk) - NkD^{t,X,p)]] 



(27r)3 ' 

which is one-electron version of more general kinetic equation for electrons in vibrating 
crystal lattice [26] . 

Thus we have demonstrated reducibility of our exact evolution equations to 
standard kinetic models. 

Unfortunately, "to break is not to make" ! The truncation of hierarchy of electron- 
many-phonon correlations is much easier than revealing of their actual predestination 
and physical sense. Therefore let us go to summing up our actual results. 

6. Discussion and resume 

The kinetic equation flSSjl takes into account only simplest mutually non-correlated 
(mutually incoherent) electron transitions, each involving one phonon. In reality, 
however, electron's interaction with any particular phonon mode k , - described by 
operators a^exp [i/cr(t)] and a[, exp [— zA;r(t)] in the Heisenberg representation over 
other modes, - depends on random electron's path r{t) determined by simultaneous 
interactions with all these modes. Hence, there are also various multi-stage electron's 
transitions which involve mutually correlated (mutually coherent) acts of absorption or 
irradiation of two, three or more phonons. Hierarchy of the above introduced correlation 
functions A„ and Dn just describes how these many-particle processes "'dress" the 
electron and constitute statistics of its random walk. 

Analogous correlated many-particle process were considered in |27] , where quantum 
amplitude of any particular electron transition was influenced by electric voltage 
fluctuations due to simultaneously realizing transitions of other electrons. There it 
was demonstrated that the correlated transitions are responsible for low-frequency 
("quasi-static", or 1/f-type) fluctuations of one-electron "transition probabilities" and 
consequently electric conductance. 

Here, we performed similar but rigorous analysis of electron (quantum particle) 
interacting with infinitely many phonon modes (boson thermostat) and showed 
possibility of shortened but complete description of this system in terms of the mentioned 
electron-many-phonon correlation functions whose arguments include electron variables 
and phonon indices only. 

Main formal result of the present paper is derivation of exact hierarchy of evolution 
equations for these correlation function and thus for electron's density matrix. That are 
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Eqs|33]or Eqsl37]and their generating functional equivalent, Eqj40l For the best of my 
knowledge, such equations never earlier were in use. 

It should be underlined that our recipe of exact shortened description can be 
extended to (generally relativistic) "dressed" electron in electromagnetic field and to 
many-electron (or, generally, many-fermion) systems. 

The correlated many-particle processes (or "kinetic events" , or "dynamic clusters" , 
etc.) in classical statistical mechanics of fluids [H El El IH El El E] and dielectric crystals 
[15] are under our attention after the work [^ where for the first time it was shown that 
a fluid particle has no definite "probability of collisions", diffusivity and mobility. In 
other words, any of such characteristics of particle's random motion undergoes scaleless 
1/f-type low-frequency fluctuations. The theory discovers them when it does not neglect 
n -particle statistical correlations which can arise between particles participating in 
chains (clusters) of n — 1 dynamically connected consecutive collisions. Infiniteness 
of variety of such chains means originality of any phase trajectory of the system [H [H] 
and impossibility of its fully adequate description in terms of two-particle events only, 
even in case of arbitrarily dilute gas [H El HI |5] . 

Close ideas concerning random motion of electrons (or holes) in electric conductors, 
- e.g. intrinsic or weakly doped semiconductors, - were formulated already in 
[ini [m EHl EH] (for later formulations see also [121 [13]), with the aim to suggest 
principal explanation of electric 1/f-noise observed in various conductors [11] and usually 
manifesting itself as 1/f fluctuations in mobilities of charge carriers [H] EOl EB 132] . 
The present work is first step from the phenomenological theory of electronic 1/f-noise 
[T0| [TT| EHl El] to its microscopic statistical-mechanical theory. 

In this respect, I would like to emphasize close conformity of the above obtained 
EqsJST] and HO] for "Brownian particle" (electron) in ideal phonon gas and previously 
obtained evolution equations for many-particle correlations "dressing" a molecular 
Brownian particle in usual classical ideal gas [H El IH E]. Solution of EqsEZl or HO] 
can be represented by an (operator-valued and branching) infinite continued fraction 
similar to that in [1]. This conformity prompts that the Eqsj37]or, equivalently, HOlalso 
hide 1/f fluctuations of diffusivity and mobility. Their visualization must be second step 
to the desired theory. 

At that, of course, none truncation of EqsEZl is appropriate. The generating- 
functional formulation, ( HOl) . of EqsJST] make it obvious that all floors of the 
corresponding continued fraction, - i.e. all the electron- ra -phonon irreducible 
correlations, - have equal rights and significance, since all they are represented in EqHOl 
by only three simple terms. And it is so at arbitrarily weak electron-phonon coupling. 
Weaker coupling implies only that upper frequency bound of 1/f fluctuations of rate of 
electron relaxation becomes lower proportionally to the rate itself, and no more. Hence, 
one needs in adequate methods of operating with the whole electron's correlation "coat" . 
This will be third step to valuable theory of electron-phonon interaction. 

I would like to acknowledge Prof. Yu. V. Medvedev for useful discussions. 
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